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GRAVITY FIELDS OF THE 
SOLAR SYSTEM 


1. INTRODUCTION 


Precise knowledge of the gravitational fields in the solar system is often an important 
ingredient in space mission planning and spacecraft design. The gravitational fields affect the 
trajectories and orbits of space vehicles and design problems involving propulsion and 
guidance systems, attitude determination and control systems, and scientific experiments. 
For example, Earth gravity models are used for calibration of launch vehicle guidance 
systems, design of scientific experiments for geodesy and general relativity, and planning for 
satellite tracking and telemetry operations. 

This monograph briefly discusses the most frequently used formulations of the gravitational 
field and defines a standard set of models for the gravity fields of the Earth, Moon, Sun, and 
other massive bodies in the solar system. These models are intended primarily for use by 
engineers and computer programmers who are not specialists in gravity field-modeling 
theory. The formulas are presented in standard forms, when possible, with instructions for 
conversion to other forms in common usage. In section 2, various formulations of the 
gravitational field are developed. The subsections are arranged in order of the sophistication 
required in their development. The first model considered is the “point source” or “inverse 
square” model, which represents the external potential of a spherically symmetrical mass 
distribution by a mathematical point mass without physical dimensions. The most obvious 
departure from symmetry in a rotating body is the latitudinal variation commonly referred 
to in terms of an equatorial bulge or a polar flattening. Accordingly, an oblate spheroid 
model is presented next. This is accompanied by an introduction to zonal harmonics. The 
spheroid model is then generalized to a representation of the field resulting from a massive 
body in terms of a spherical harmonic expansion. The latter formulation is the basis for a 
number of the spherical harmonic models which have been developed for the Earth and 
Moon. These models and their application to NASA missions are discussed. In addition, the 
triaxial ellipsoid model is presented because of its suitability for modeling lunar gravity and 
use in lunar missions. 

Section 3 provides guidance in selection of gravitational models for the Earth and Moon and 
gives values of the basic parameters used to describe the gravity fields of the other bodies in 
the solar system. 

Appendix A defines the symbols used in this monograph. 

2. STATE OF THE ART 

The gravitational field in the vicinity of a celestial body may be described in several ways. In 
general, planets and other massive celestial bodies are not perfectly symmetrical in shape 
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and their mass is not uniformly distributed. Consequently, in each case complete 
representation of the gravitational field would require an infinite set of orthogonal functions 
such as spherical harmonics. If the shape and mass distribution are reasonably regular, 
however, the field may be approximated by much simpler representations. A common 
simplifying technique is to assign the celestial body a regular geometrical shape that closely 
approximates the true shape and to assume a uniform density for the material of the body. 
This approach results in a gravitational field model whose equipotential surfaces are shaped 
approximately like the body’s surface. The usual shapes assumed are a sphere, an ellipsoid of 
revolution, or a triaxial ellipsoid. A spherical body has spherically-symmetrical equipotential 
surfaces. The equipotential surfaces external to an ellipsoid of revolution possess axial 
symmetry but display a latitudinal (zonal) variation in radius similar to that of the body’s 
surface. The equipotential surfaces external to a triaxial ellipsoid vary in radius with 
longitude and latitude in a manner similar to the meridional and equatorial ellipticities of 
the body’s shape. With the assumption of uniform density, a body of arbitrary shape results 
in equipotential surfaces which can be represented by an infinite sequence of spherical 
harmonics. 

The gravitational potentials that result from the foregoing models are presented in the 
following subsections. 

2.1 Point-Mass Representation 

The simplest way to express the external gravitational field of a celestial body is to treat it 
as a spherical, uniform body or one whose mass can be considered concentrated at a point. 
With this assumption, the field is described by a single parameter, the gravitational constant 
M, which is the product of the universal gravitational constant G and the mass of the body 
M. For all the major bodies in the solar system except Earth, Moon, and Mars, the only 
parameters of the gravitational fields which may be considered well-known are the p values 
and the ratios of the mass of the Sun to the masses of the bodies. The mathematical 
premises for p and M are given in appendix B. 

2.1.1 Application of Point-Mass Model 

On the basis of Newton’s universal law of gravitation, the point mass representation of a 
celestial body results in a potential field that is spherically symmetric about the point mass. 
The location of the point mass is the center of mass of the body, the only point in the body 
through which the application of any external force will result only in translational motion. 
The attractive force F that results from such a model is given by the inverse square law. The 
force on a unit mass located a distance R from the center of mass of the body is given by 


This force may also be expressed as the magnitude of the gradient of a scalar gravitational 
potential U by 



( 2 ) 
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If the potential of a celestial body is needed at a distance R that is large compared to the 
body’s diameter or if precise knowledge of the effect of small perturbing forces is not 
required, the potential may be adequately represented with the point-mass model by 
specifying the values of ju and R. 

For a spacecraft traversing cislunar space, the gravitational potential fields of the Earth and 
Moon may be adequately described with a sum of point-mass models for the Earth and 
Moon by specifying the value of ju for each body and the distances R e and R m from the 
Earth and Moon, respectively. In interplanetary space, far from any massive body, the 
point-mass models for the major planets, the asteroids, and the Sun give a sufficient 
description of the potential. In the case of a planet with natural satellites, a good 
approximation of its gravitational field at large distances can be obtained by considering the 
planet and satellites as a single mass point located at the center of their combined mass with 
a mass equal to their combined mass. In such cases, the consideration of the natural 
satellites is particularly important in generating accurate planetary ephemerides. When 
ephemerides are generated over long time periods, even such small corrections are 
important. 

For most of the planets and all the minor celestial bodies, the point-mass model is the only 
one that can presently be used with confidence. In some cases, notably Pluto, even the ratio 
of its mass to that of the Sun and the ju value are still quite uncertain. Knowledge of the 
detailed effects of the gravity fields of Mars and Venus on orbiting vehicles is becoming 
important, however, with the advent of orbiting and soft-landing missions such as Viking. 


2.1.2 Gravitational Constant and Mass 


2.1. 2.1 Sun 


Table 1 gives recent determinations of the gravitational constant of the Sun, The values 
given by references 1 and 2 are extremely close to each other, and both are generally 
accepted in gravitational work. The Astronomical Ephemeris value is presented for historical 
perspective. It was obtained without benefit of data from cislunar and interplanetary 
spacecraft missions. Because the uncertainty in the Jet Propulsion Laboratory (JPL) value 
(ref. 2) is small, this value is adopted herein. 


The value for the Sun’s mass, M s that corresponds to the foregoing adopted value for n s is 
obtained by 
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TABLE 1 

GRAVITATIONAL CONSTANT FOR THE SUN 


Source ! 

M s (km 3 /s 2 ) 


Astronomical Ephemeris (ref. 3) 
Anderson et al (ref. 1 ) 

JPL Value (ref. 2) 

1.3246 x 10 11 

1.327125 x 10 11 

1.32712499 x 10 11 ± 15x 10 3 


where G is the universal gravitational constant, equal to 6.668 ± .005 X 10' 23 km 3 /g s 2 . 
The resulting value of 1.9903 X 10 33 g is adopted for M s . 


2. 1.2. 2 Planets 

The usual method of presenting the mass of a planet is in terms of the reciprocal mass ratio, 
(Mp)' 1 where M^ is the observable quantity M p /M s (appendix B). Tables 2 through 9 give 
the best estimates presently available for the reciprocal mass ratios for each of the planets 
except Earth. The form of the presentation has been adopted from reference 4. The U.S. 
Naval Observatory (USNO) values are weighted means of all the recent determinations done 
at the Observatory and elsewhere. The JPL values were adopted in 1968 and taken from 
reference 2. The Massachusetts Institute of Technology (MIT) values are based on a larger 
data set than any of the others and include optical and radar observations. Kovalevsky’s 
approach was similar to that of USNO although a smaller set of determinations was used. 
Kovalevsky’s values are presented in terms of a range of equally probable values. Anderson’s 
results are from a recent analysis of Mariner 5 data. 

The reciprocal mass ratios for the major planets are given in tables 2 through 9. The 
uncertainty values therein that were taken from references 2, 4, and 5 are computed from 
formal statistics. They show the dispersion of data but not necessarily the probable error in 
derivation of planetary masses. 

On the other hand, the uncertainty values accompanying references 6 through 16 in tables 2 
through 9 are estimates of the real error in the planetary mass determinations. Accordingly, 
the real error estimates serve as the basis for the weighted mean values given in the tables for 
the planetary reciprocal mass ratios. 
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TABLE 2 

RECIPROCAL MASS RATIO OF MERCURY 


Source 

(M p ) 1 

Uncertainty 

Weighted Mean Value 
USNO (ref. 4) 

JPL (ref. 2) 

MIT (ref. 6) 
Kovalevsky (ref. 7) 
Ash (ref. 8) 

Howard (ref. 9) 

6,021,900 

5.987.000 

5.983.000 

6.020.000 

5.900.000 t 

6.025.000 
6,023,600 

32.000 

25.000 

10.000 
o 6,100,000 

15,000 

600 


TABLE 3 

RECIPROCAL MASS RATIO OF VENUS 


Source 

(M' p ) _1 

Uncertainty 

Weighted Mean Value 

408,522.7 


USNO (ref. 4) 

408,519 

11 

JPL (ref. 2) 

408,522 

3 

MIT (ref. 6) 

408,522 

3 

Kovalevsky (ref. 7) 

408,512 to 408,532 

Mariners 2 and 4 (ref. 10) 

408,521.8 

1 

Anderson (ref. 1 1 ) 

408,523.5 

1 

Howard (ref. 12) 

408,523.9 

1.2 


TABLE 4 

RECIPROCAL MASS RATIO OF MARS 


Source 

(M'p)- 1 

Uncertainty 

Weighted Mean Value 
USNO (ref. 4) 

3,098,710 

3,098,709 

9 

JPL (ref. 2) 

3,098,700 

100 

MIT (ref. 6) 

3,098,700 

30 

Kovalevsky (ref. 7) 

3,098,650 to 3,098,750 

Mariner 4 (ref. 13) 

3,098,708 

9 

Mariner 9 (ref. 14) 

3,098,720 

70 


TABLE 5 

RECIPROCAL MASS RATIO OF JUPITER 




Source 

(M'p)' 1 

Uncertainty 

Weighted Mean Value 
USNO (ref. 4) 

JPL (ref. 2) 

MIT (ref. 6) 

Kovalevsky (ref. 7) 

Jupiter Monograph (ref. 15) 
Pioneer 10 (ref. 16) 

1,047,373.6 

1,047.366 

1,047.3908 

1,047.4 

1,047.34 t( 

1,047.39 

1,047.342 

0.007 

0.0074 

0.1 

1,047.39 

0.04 

0.02 
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TABLE 6 

RECIPROCAL MASS RATIO OF SATURN 


Source 

(M p ) 1 

Uncertainty 

Weighted Mean Value 
USNO (ref. 4) 

JPL (ref. 2) 

MIT (ref. 6) 
Kovalevsky (ref. 7) 

3,498.5 

3.498.1 

3.499.2 
3,498.5 

3,497 t 

0.4 

0.4 

0.5 

o 3,500 


TABLE 7 

RECIPROCAL MASS RATIO OF URANUS 


Source 

(Mp ) 1 

Uncertainty 

Weighted Mean Value 
USNO (ref. 4) 

JPL (ref. 2) 

MIT (ref. 6) 
Kovalevsky (ref. 7) 

22,920 

22,800 

22,930 

22,900 

22,600 t 

107 

6 

200 

0 23,000 

1 


TABLE 8 

RECIPROCAL MASS RATIO OF NEPTUNE 


Source 

(M p ) 1 

Uncertainty 

Weighted Mean Value 
USNO (ref. 4) 

JPL (ref. 2) 

MIT (ref. 6) 
Kovalevsky (ref. 7) 

19,323 

19,325 

19,260 

19,400 

19,200 t 

1 

26 

100 

100 

19,400 

1 


TABLE 9 

RECIPROCAL MASS RATIO OF PLUTO 


Source 

(Mp ) 1 

Uncertainty 

Weighted Mean Value 
USNO (ref. 5) 

JPL (ref. 2) 

MIT (ref. 6) 
Kovalevsky (ref. 7) 

1.934.000 
3,000,000 

1.812.000 

3.500.000 

1 .500.000 t 

500,000 
40,000 
2,000,000 
o 2,500,000 
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The gravitational constant of a planet, ju p may be computed from 

«p=gm,m;= ( .,/(m;)- 1 o) 

where (M p ) _1 is found in tables 2 through 9. The advantage of computing p p in this way is 
that it avoids the uncertainty in our knowledge of G. 


2. 1.2. 3 Earth and Moon 

The usual means of expressing the mass of the Earth is in terms of the inverse ratio of the 
combined mass of the Earth and Moon to the mass of the Sun: 




'M + M ' 

e m 


M 


-1 


(4) 


Table 10 contains the best available estimates of (M* V 1 . The comment in section 2. 1.2. 2 
concerning the uncertainty ranges of reciprocal mass ratios of planets for different 
investigators applies also to values from these investigators in table 10. 


Table 1 1 contains values of p e = GM e that are computed by three different methods. With a 
spheroid model for the Earth (sec. 2.2), ju e can be estimated in terms of the Earth’s rotation 
rate and the dynamic flattening. A second method uses radar values for the mean distance of 
the Moon in a modified form of Kepler’s equations. The third determination is from the 
tracking of lunar probes by JPL. Also included in table 11 are a value recently determined 
by Esposito and Wong from Mariner 9 data (ref. 17) and the values adopted by JPL in 1968 
(ref. 2). 

The masses of the Earth and Moon can then be separated using the ratio M e /M m , given in 
table 12. 


Prior to the launching of spacecraft into cislunar and interplanetary space, the ratio of the 
mass of the Moon to that of the Earth was known to three significant figures, 
M e /M m =81.3 (ref. 21). This value may still be used when only three significant figures are 
required. The Earth-Moon system parameters that can be estimated directly are the 
gravitational constant u = GM and the mass ratio M /M m . As in the case of the mass of 
the Sun (sec. 2. 1.2.1), the absolute value of M m is known only as accurately as G, and G is 
known only to about three significant figures (sec. 3.1). 

Recent determinations of u and M / M are listed in table 12. The values for Rangers 6, 7, 
8, and 9 and Mariners 2 and 4 were computed from reference 21; Blackshear’s value 
represents a best fit to his latest gravity model derived from Lunar Orbiter data; and the 
values derived from Mariners 5, 6, and 7 and Pioneers 8 and 9 were communicated by Null 
of JPL and represent JPL’s adopted values as of 1969. 
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TABLE 10 

RECIPROCAL MASS RATIOS FOR EARTH PLUS MOON 


Source 

(M' I' 1 

e-m 

Uncertainty 

Weighted Mean Value 
USNO (ref. 4) 

JPL (ref. 2) 

MIT (ref. 6) 
Kovalevsky (ref. 7) 

328,900.21 

328,900.12 

328,900.1 

328,900 

328,900 t 

: 

0.20 

0.3 

1.0 

0 328,930 

1 


TABLE 11 

THE EARTH'S GRAVITATIONAL CONSTANT 


Source 

M e (km 3 /s 2 ) 

Dynamic Flattening (ref. 18) 

Mean Distance of the Moon (ref. 19) 
Lunar Probes (ref. 20) 

Mariner 9 (ref. 17) 

JPL Adopted Value (ref. 2) 

398,603.2 

398,600.1 

398,600.9 

398,600.8 

398,601.2* 


* Adopted for this monograph. 


TABLE 12 

GRAVITATIONAL CONSTANT AND MASS RATIO OF THE MOON 


Source 

At m (km 3 /s 2 ) 

M /M 

e m 

Weighted Mean Value 

4902.78 

81.30090 

Rangers 6, 7, 8, and 9 (ref. 21) 

4902.65 ±0.16 

81.30245 ± .00246 

Mariners 2 and 4 (ref. 21) 

4902.735 ± 0.21 

81.30175± .00315 

Blackshear 1 

4902.867 

81.2994 

Mariners 5, 6, and 7 , f 7 o\ 
Pioneers 8 and 9 1 ' 

4902.801 ± .022 

81.30071 ± .00036 

JPL Adopted Value (ref. 2) 

4902.78 +0.06 

81.3010 +.001 


1 Personal Communication, August 28, 1970. 
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2.2 Spheroid Model 

2.2.1 Flattening and Zonal Harmonics, J 2 and J 4 

A rotating celestial body may be approximated with good accuracy by an ellipsoid 
(spheroid) of revolution. The gravitational field of such a body has axial symmetry but 
varies with latitude because of its ellipsoidal shape. In most formulations, this variation is 
expressed in terms of the second and fourth even-zonal harmonics of the spherical harmonic 
expansion having constant coefficients, J 2 and J 4 . As shown in appendix C, these 
coefficients may in turn be expressed in terms of the rotation and the dynamical flattening f 
which is given by 


f = 


( a e - b ) 


a 


e 


where a e = the equatorial radius 
b = the polar radius 

The departures from the spherically symmetric field that result from the J 2 and J 4 terms 
may be represented geometrically by relative highs and lows in the equipotential surfaces as 
compared to a sphere (fig. 1 ). 

The geometrical flattening of the major planets has been estimated from optical 
observations of their shapes. From these estimates and the expressions relating the rotation, 
J 2 , and J 4 to f, J 2 and J 4 may be estimated. For planets with natural satellites (particularly 
Mars, Jupiter, and Saturn), optical tracking of the satellites can be used to estimate the 
magnitude of the secular perturbations from which J 2 and J 4 may be computed. 

The gravitational field that results from a spheroidal body is treated in references 23, 24, 25, 
and 26; references 23 and 27 give a brief mathematical description of the gravitational field 
of a spheroid body with a discussion of Legendre functions. 


2.2.2 Application 

The spheroid model is an important tool in formulating the gravitational field of a planet. 
The parameters that describe it can be derived directly from observations of orbiting 
satellites which show large secular changes in the right ascension of the ascending node £2 
and in the argument of periapsis co. Because of the large secular motions of close-Earth 
satellites, the spheroid model is the simplest model that can be recommended for orbit 
determination; however, more detailed and accurate models are usually used. The ellipsoidal 
equipotential surface of the spheroid model also serves as the common reference surface for 
many geodetic applications. 
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2.2.3 Major Planets 


Data from tracking of natural satellites and spacecraft in the vicinity of the nearer planets 
have yielded estimates of only J 2 and J 4 for some of them. In some cases, the dynamical 
oblateness can be estimated only from the optically-observed geometrical flattening (sec. 
2.2.1). Such estimates, however, require assumptions about the plasticity, viscosity, and 
modulus of elasticity of the planetary interiors that are often little more than educated 
guesses. 




+ High regions 
- Low regions 


Figure 1.— Highs and Lows of a Spheroidal Equipotential Surface 
Relative to a Spherical Surface 



2.2.3. 1 Mercury and Venus 

Mercury and Venus have no significant flattening (ref. 3). As Venus has no natural satellites, 
the only available source of data on Venus’ dynamic oblateness is the Mariner 5 flyby from 
which Anderson (ref. 1 1) has estimated 

J 2 (Venus) = 2.7 ± 0.9 X 1 O' 5 

No estimates have been made for J 4 . No dynamical estimates of J 2 and J 4 have been made 
for Mercury. 

2.2.3.2 Mars 

J 2 for Mars has been estimated from the motion of its two natural satellites by Wilkins (ref. 
28) and from Doppler tracking of Mariner 4 by Null (ref. 13). Because the actual quantity 
derived in both cases was not J 2 but the product J 2 where R m is an assumed value for 
the equatorial radius of Mars, the computed value of J 2 is affected by the determination of 
R . The values for J„ in table 13 (refs. 13 and 28) assume a value for R„ of 3394 km in 
contrast to the radius of 3402 ± 8 km for this monograph (sec. 3.2.2). Wilkins result has 
been adjusted to reflect a value of 3394 km for R m instead of his original value of 3409 km. 


TABLE 13 
J 2 FOR MARS 


Source 


Natural Satellites (ref. 28) 
Mariner 4 (ref. 13) 

Lorell & Shapiro (ref. 14) 
Jordan & Lorell (ref. 29) 

.001968 ± .000006 
.00187 +.00007 
.00196 ±.00001 
.001964 ± .000006 


It was not feasible to estimate J 2 from the motion of Mariners 6 and 7 because of 
nongravitational perturbations such as gas venting that acted on the spacecraft. J 4 has not 
been estimated for Mars. The flattening for Mars has been estimated geometrically to be 
(192)' 1 (ref. 3); this compares with Wilkins’ estimate of (190.4 ± 1.9)" 1 and estimate of 
(190.8 ± .7)' 1 by Lorell and Shapiro for the dynamical flattening (ref. 28). 

2.2. 3.3 Jupiter 

A value for the geometrical and dynamical flattening of Jupiter are computed in reference 
13 as (16.35)' 1 . Anderson et al. (ref. 16) calculate dynamical flattening as (15.456 ± 0.24) 
from Pioneer 10 results and give best available determinations of J 2 and J 4 : 

J 2 (Jupiter) = (147.2 ±0.4) X 10' 4 
J 4 (Jupiter) = -(6.5 ±3.8) X 10' 4 
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2. 2.3.4 Saturn 


In reference 30, the geometrical and dynamical flattenings of Saturn are computed to be 
(9.5 ± 0.7)' 1 and (10.3 ± 0.5)' 1 , respectively. The values of J 2 and J 4 computed in 
reference 26 are 

J 2 (Saturn) = (166.5 ± 2.0) X 10' 4 
J 4 (Saturn) = (9.6 ± 1.0) X 10' 4 

Reference 31 notes that these values contain a contribution of undetermined magnitude 
from Saturn’s rings. 

2.2. 3.5 Uranus, Neptune, and Pluto 

Reference 32 contains estimates for the geometrical flattening of Uranus and Neptune of 
(33 ± 33)' 1 and (50 ± 50)' 1 , respectively. This notation indicates that the correct values 
may lie anywhere in the range of expected uncertainties with virtually equal probability. 
Reference 33 gives estimates of the geometrical flattening of (18)' 1 for Uranus and (60)' 1 
for Neptune, and a value of J 2 for Neptune of approximately 0.005. 

No estimates of Pluto’s spheroid parameters have been made. 

2.2.4 Sun 

Determination of the Sun’s zonal harmonics presents a different problem from that 
encountered with the planets. Because of the large mass of the Sun, the effects of general 
relativity must be accounted for in determining J 2 and J 4 from such observed quantities as 
the advance of the perihelion of Mercury. The observed rate of advance agrees well with the 
prediction of general relativity but leaves the contribution from the Sun’s oblateness quite 
uncertain. Anderson (ref. 1) has estimated that J 2 for the Sun has an upper limit of 10' 5 
with an uncertainty of 70 percent. Dicke (ref. 34) has postulated a certain amount of 
geometrical flattening of the Sun in an attempt to support his own relativistic theory that 
would require that the Sun’s oblateness contribute about 10 percent of the motion of 
Mercury’s perihelion. His value for f was 5.0 ± 0.7 X 10' 5 . A recent determination by 
Oesterwinter (ref. 35) did not estimate J 2 precisely for the Sun but gave it an upper limit of 

1.4 X 10' 5 . 

2.3 Spherical Harmonic Model 

The representation of a planet’s surface topography in spherical coordinates is appropriate 
because most planets are nearly spherical. The formulation of the planet’s gravitational 
potential in terms of spherical harmonic expansion follows naturally because the 
gravitational field is strongly related to the planet’s shape. In recent years spherical 
harmonic models have been constructed for the gravity fields of the Earth and Moon; and as 
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data becomes available, attempts will be made to construct models for other planets. 
Appendix D gives the theoretical formulation of the gravitational potential in terms of 
spherical harmonics and discusses the physical significance of their mutual orthogonality. 


2.3.1 Orthogonal Functions 

All irregularities in shape and mass distribution contribute to the total gravitational 
potential of a body. To model such irregularities precisely, a formulation is required which 
uses an infinite sequence of orthogonal functions. The sequence of orthogonal spherical 
harmonics is appropriate for representation of a function over a spherical surface which is 
analogous to the use of Fourier series for functions in a rectilinear space. 


2.3.2 Representation of the Gravitational Potential 

2.3.2. 1 Zonal, Sectorial, and Tesseral Harmonics 

Spherical harmonics may be visualized as small, periodic adjustments to a perfectly spherical 
shape that result in a surface which oscillates above and below the spherical surface at 
regular intervals; positive values represent local “highs” relative to the reference sphere and 
negative values, “lows.” The categorization of the harmonics is based on the geometrical 
pattern of the highs and lows. Zonal harmonics possess longitudinal symmetry and oscillate 
in sign only with latitude. Sectorial harmonics possess latitudinal symmetry and oscillate in 
sign only with longitude. Tesseral harmonics comprise all the other harmonics of varying 
degree and order and change in sign over a latitude/longitude grid of “tesserae.” An 
analytical description of the foregoing concepts is given in appendix E. 

2.3.2.2 Resonant Orbits 

Certain Earth satellite orbits, described as “resonant,” experience large perturbing effects 
from particular harmonic coefficients (ref. 26). A resonant orbit has a mean motion 
commensurate with Earth’s rotation, i.e., the ratio 24/P is a rational number where P is 
approximately the satellite period in hours. The most pronounced resonances occur when 
24/P is an integer. During each commensurate period, the geographic trace of the satellite 
repeats, so any anomaly in the geopotential, however small, has an opportunity to build up 
small perturbations into large displacements as long as the commensurability is maintained. 
Such orbits have been recognized since 1960 (ref. 36) as the reasonable ones to observe for 
the determination of certain of the longitude-dependent terms of the geopotential. 
Appendix E gives a more detailed discussion and a list of harmonic coefficients so 
determined. 


2.3. 2. 3 Normalization and Conversion 

The coefficients C nm and S nm which appear in the expression for the geopotential in 
spherical harmonics, often are normalized or combined with other coefficients to yield 
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coefficients that have clear physical interpretations. Appendix F gives the generally adopted 
normalization procedure and relates the various forms of geopotential expression to the 
form adopted herein (equation E-2). 


2.3.3 Current Models 

2.3.3. 1 Earth 

Many spherical harmonic models have been developed for Earth gravity since the first 
artificial Earth satellite flew in 1959. Each successive satellite launch and each advance in 
tracking system accuracy have led to more complete and accurate models representing a 
wider variety of orbits. Notable models that have led to present state of the art models are 

• The 4 by 4 and 6 by 6 models of Izsak derived from optical data (refs. 37 and 38) 

• The 4 by 4 and 6 by 6 Doppler data models of Guier (ref. 39) and Anderle (ref. 40) 

• The 8 by 8 model of Kaula (ref. 41) based on both optical and surface gravity data 

• The 8 by 8 Doppler data model of Guier and Newton (ref. 42) 

• The SAO-M 1 model based on optical and surface gravity data described by Gaposchkin 
(ref. 43) and Lundquist (ref. 44) 

• The SAO-69 model based on optical and surface gravity data (ref. 45) 

• The GEM-6 gravity model (ref. 46) 

The Goddard Earth Model-6 (GEM-6) (ref. 46) reflects a denser distribution of tracking 
stations, more different types of tracking data, and a wider range of orbit inclinations than 
any other such model and includes observations of 26 satellites and considerable amounts of 
surface gravity data. The GEM-6 model consists of a complete set of tesseral and sectorial 
harmonics through degree and order 16, zonal harmonics through degree 22, and other 
selected resonant terms (sec. 2.3. 2. 2) of orders 9, 11, 12, 13, and 14. In deriving the 
coefficients of the GEM-6 model, both optical and laser tracking data were used. GEM-6 is 
recommended for the general applications described in section 3.3.2. 1. 


2. 3. 3.2 Moon 

The spherical harmonic expansions of the gravitational potential of the Moon are 
incomplete and less reliable than those for Earth because of a lack of tracking data from a 
broad distribution of orbit inclinations and lack of information from behind the Moon. 
However, interesting results have come from studies of Apollo flights and Explorer 49 (refs. 
47 through 50). One spherical harmonic expansion, the LI model of NASA Langley 
Research Center, was used for the Apollo 11 through 17 missions (ref. 51). Although it is 
difficult to measure quantitatively the accuracy of global models derived from the Lunar 
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Orbiter data, it may be concluded that they model the data from which they were derived 
very well. Therefore, the models could be considered in analysis of satelite orbits similar to 
those flown by the Lunar Orbiters. Data from the Apollo flights have generally not been 
used for spherical harmonic models because Apollo tracking data are corrupted by effects 
such as gas venting. 

No spherical harmonic expansion models are currently recommended for all Moon missions. 
The triaxial ellipsoid model (sec. 2.4) is better suited for a wide variety of lunar applications 
except for low altitude lunar orbiters. 

2.3.4 Model Truncation 

This section provides guidance for the application of a spherical harmonic model to 
particular requirements. The spherical harmonic models that have been developed for the 
Earth’s gravitational potential contain a finite number of terms. Such models are usually 
most effective for orbit prediction when applied to orbits similar to those used in deriving 
the models. 

Because of computational limitations, users of gravity models often require models 
containing only a small number of terms. Two options may be possible; a model of the 
desired size may already be available, but more often the user will have to reduce a larger 
model. In the latter case, the problem of which terms to retain still must be met. 

If a user requires a smaller model than the GEM-6 model, the following guidelines may assist 
in selecting which terms to retain: 

• All terms which are at or near resonance for the application orbits 

• Only low-degree terms at higher altitudes because of the l/r n dependence in the 
potential where n equals the degree of the term 

More detailed quantitative rules for determining which terms are to be retained or how 
much accuracy is lost in truncation are given in appendix G. 

2.3.5 Model Accuracy and Gravimetry 

Model accuracy in the context of space applications refers to the ability of a model to 
predict the position of a satellite in a known orbit. A prediction accuracy of ± 30 meters is 
widely quoted as being generally attainable. 

Such accuracy designations, however, are misleading and often lead to misconceptions. The 
accuracy with which the motion of a satellite can be predicted depends on the length of 
time over which the prediction is made, parameters of the initial orbit, the method of 
computation, the way in which the results are expressed, and the manner in which the 
model is applied. It is, therefore, impossible to assign a single number to represent the 
accuracy of a model. 
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For the GEM-6 gravity model, reference 45 provides a satisfactory discussion of model 
accuracy. 

Surface gravimetry has recently been investigated as a possible standard for evaluating the 
accuracy of satellite models. It has been found that over small areas, gravimetric data can 
yield far more accurate gravity values than can be deduced from satellite data. Thus a 
well-measured portion of the geoid can serve as a reference standard for satellite models. 
The severe limitation on surface gravimetry is that neither the technological nor the 
economic means have been available to obtain data over much of the ocean surface. 

The use of satellite data involves the reverse situation. Although satellite data represent the 
integrated effects of gravity from every point on the Earth’s surface, these measurements 
lack fine resolution. A model complete to degree 30 has a resolution element of 180/30 or 6 
degrees of arc relative to the Earth’s center. When applied to surface gravity, such a model 
can yield average values of gravity over grid elements no smaller than 360 nautical miles (6 
degrees of arc) on a side. Where surface gravity data are available, however, it is common to 
attain accurate average gravity over surface regions extending 1 degree of arc (60 nautical 
miles) or smaller. 

Thus, one of the major problems encountered in using surface gravity is the selection of an 
appropriate grid size over which comparisons with satellite gravity can be made. A surface 
gravity grid contains too many details to be approximated by satellite gravity. Moreover, 
because of the limited coverage of gravimetry data, it cannot be extrapolated realistically 
over the entire satellite resolution element. Kaula (ref. 52) selected a 300-nautical mile (5 
degrees) grid size over which to compare surface measurements with data from some early 
satellite models. 

2.4 The Triaxial Ellipsoid Model (Moon) 

2.4.1 Principal Moments of Inertia 

The triaxial ellipsoid is a refinement of the spheroid of revolution model (sec. 2.2). The 
spheroid of revolution yielded a gravitational field which possessed longitudinal symmetry 
but reflected a latitudinal ellipticity. In physical terms, the principal moments of inertia of 
the spheroid of revolution, A, B, and C, defined with C measured about the rotation axis, 
satisfy A = B C. The spheroid model potential, as formulated in the general case trom 
MacCullagh’s formula (ref. 20), is 


p A + B + C-3I 
U =— + G + O 

r 2r 3 


( 7 ) 


where r - the distance from the center of mass to an external point P 

I = the moment of inertia about the line joining the center of mass and P 
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If the spheroid potential (sec. 2.2) is compared to this expression for U, J 2 can be written in 
terms of the principal moments of inertia by 


“a l 


where a g = the equatorial radius. 

For the triaxial ellipsoid where A ¥= B, the corresponding expression for J 2 is given by 

T = c 4 (a + B) 

J 2 ^ 

Ma 2 

e 


Where a e = the mean equatorial radius. 

In this model there is an equatorial as well as a polar flattening. As the polar flattening is 
given principally by the spherical harmonic coefficient C 20 (equals -J 2 ), the equatorial 
flattening is represented by the coefficient C 22 , which is related to the moments of inertia 
(refs. 26 and 53) by 



Derivations pertaining to the preceding discussion of the potential from a triaxial ellipsoid 
may be found in references 24 and 54. The gravity potential of a triaxial ellipsoid may also 
be expressed in terms of the spherical harmonic expansion by using only three terms, the 
central force term (ju), C 20 , and C 22 - 

2.4.2 Application 

The triaxial model has been applied to the Earth, but it has not yielded significant 
improvements over the spheroid model. The main difference between these two models is 
the C 22 term. Many investigators, including Wagner (ref. 55), have studied the ellipticity of 
the Earth’s equator and found that C 22 is nearly three orders of magnitude smaller than 
C 20 and is therefore negligible for applications requiring less detail than the spherical 
harmonic expansion models presented in section 2.3.3. 1. For the Earth, this result can be 
qualitatively expressed in terms of principal moments of inertia by 

| C - A | » | C - B I > I A-B | 

The triaxial ellipsoid is more useful for the Moon. From the moments of inertia ratios, given 
in reference 56, it can be inferred that C > B > A for the Moon and further that 
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and 


C-A 


B- A 


= 2.79 ± .04 


C-B 


B- A 


1.79 ± .04 


Both of the foregoing ratios are about 500 for the Earth, i.e., the equatorial flattening is 
negligible compared to the polar flattening. Because of the more significant equatorial 
flattening on the Moon, the triaxial ellipsoid is an appropriate model and often has been 
used to describe the lunar gravitational field. Recent determinations of C 20 and C 22 for the 
Moon are given in table 14. 


TABLE 14 

TRIAXIAL COEFFICIENTS FOR THE MOON 


Source 

^2 0 

c 22 

Melbourne et al. (ref. 2) 
JSC 1 

Bender et al. (ref. 56) 

! 

-2.0711 xIO 4 
-2.07108 x 10‘ 4 
-2.04 xIO 4 

0.20716 x 10 4 
0.20716 x 10 4 
0.223 x 10' 4 


1 W. Wollenhaupt, JSC, Personal Communications, August 1970. 


3. CRITERIA 


The descriptive parameters given in this section should be used to establish reference gravity 
fields for space mission planning and the design of space vehicles, experiments, and 
instrumentation. 


3.1 Sun and Astronomical Unit 

The value of the Sun’s gravitational constant ju recommended for use is 

p s = 1.327125 X 10 1 1 km 3 /s 2 
The corresponding value of the Sun’s mass M s is 

M s = 1.9903 X 10 33 g 
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which results from 


M s = 




G 


where G is the universal gravitational constant, taken as 6.668 ± .005 X 1(T 23 km 3 / g s 2 . 

The parameter is the only well-known gravitational parameter for the Sun. Of all the 
other parameters, only J 2 has been estimated, and then only an upper limit of 10' 5 is 
recommended. The same limit may be applied to the geometrical flattening f. 

The value recommended for the Astronomical Unit (AU) is 149,597,893 km. The 
International Astronomical Union (IAU) value of 149,500,000 km is important for 
interpretation of historical astronomical observations. 


3.2 Planets 

3.2.1 Masses and Gravitational Constants 

The recommended ratios of the mass of the Sun to those of the planets, (Mp" 1 , and the 
gravitational constants of the major planets p p are given in table 15. 

3.2.2 Planetary Radii and Mean Distances From Sun 

Table 16 gives the radii of the major planets and the mean distances from the Sun to each 
planet. 


3.2.3 Harmonic Coefficients 

At present, the only harmonic coefficients that have been reliably estimated for the major 
planets are the spheroidal harmonics J 2 and J 4 which are used in the model presented in 
section 2.2.2. In some cases, all that is available is an estimate of the observed geometrical 
flattening. Table 17 gives the recommended values for the parameters that have been 
estimated. 
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TABLE 15 

RECIPROCAL MASS RATIOS AND GRAVITATIONAL 
CONSTANTS FOR THE PLANETS 


Planet 

M p (km 3 /s 2 ) 

(M'p)' 1 

Mercury 

Venus 

Mars 

Jupiter 

Saturn 

Uranus 

Neptune 

Pluto 

22,032* 

324,860 

42,828 

126,709,801 

37,934,115 

5,790,249 

6,868,111 

68,621 

i 

6,021,900 

408,522.7 

3,098,707 

1,047.3736 

3,498.5 

22,920 

19,323 

1,934,000 


* Reference 9. 


TABLE 16 

RADII AND MEAN DISTANCES 
FROM THE SUN OF THE MAJOR PLANETS 


Planet 

Equatorial 

Radius 

(km) 

Mean Distance From Sun 1 

(AU) 

(km) 

Mercury (ref. 9) 

2,439 ± 1 

.387099 

57,909,195 

Venus (ref. 57) 

6,052 ± 6 

.723332 

108,208,943 

Earth (refs. 2 and 58) 

6,378.16 ± .005 

1.0 

149,597,893 

Mars (ref. 57) 

3,402 ± 8 

1.523691 

227,940,963 

Jupiter (ref. 27) 

71,422 ±200 

5.202803 

778,328,366 

Saturn (ref. 30) 

59,800 ± 350 

9.538843 

1,426,990,814 

Uranus (ref. 32) 

27,000 ± 1,000 

19.181951 

2,869,579,453 

Neptune (ref. 32) 

25,200 ± 200 

30.057779 

4,496,580,407 

Pluto (ref. 57) 

2,250 ± 1,150 

39.43871 

5,899,947,919 


1 Reference 3. 


TABLE 17 

SPHEROID MODEL COEFFICIENTS FOR THE PLANETS 


Planet 

Geometrical 
Flattening, f 

CN 

— D 

J 4 

Mercury 

0.0 



Venus 

0.0 

~ 8.4 x 10' 6 


Mars 

( 1 92)' 1 

1.9 x 10 3 


Jupiter 

( 1 5.456)' 1 

147.2 x 10' 4 

-6.5 x 10' 4 

Saturn 

(9.5)' 1 

0.017 

9.6 x 10' 4 

Uranus 

(25) 1 



Neptune 

(50) 

~ 0.005 
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3.3 EARTH 


3.3.1 Mass and Gravitational Constant 

In planning for deep space trajectories (far from the Earth and Moon), the Earth and Moon 
should be treated as a single body whose mass M is given by 

M s 

M e-m 328,900 

where M s is the mass of the Sun. The gravitational constant of the Earth-Moon p e m is 
403,504 km 3 /s 2 . 

In or near cislunar space the Earth and Moon should be treated as separate bodies. The mass 
and gravitational constant of the Earth are given by 

M s 

M e 332,945.4 
M e =398,601.2 km 3 /s 2 

3.3.2 Spherical Harmonic Models 

The gravitational field of a rotating celestial body is generally that of a sphere modified by 
effects of rotation and irregularities in mass density and topography. Such a field can be 
described effectively by a spherical harmonic expansion model. 


3. 3. 2.1 General Application 

The GEM-6 model of the Earth’s gravitational potential should be used when there is a 
requirement involving many different satellite orbits or an undefined or unspecified satellite 
orbit. 

Table 18 gives the values of the coefficients of the GEM-6 model. Because of its size, the 
GEM-6 model may be unwieldy or inefficient for certain applications. In such cases, its size 
may be reduced in accordance with the guidelines presented in section 2.3.4. 


3.3.2. 2 Orbits at Resonant Altitudes 

When a satellite is to be flown at or near one of the resonant altitudes (sec. E.4, appendix 
E), the mission planner should select the associated resonant coefficients from an 
appropriate model and use them in the GEM-6 model. For example, if the period of an orbit 
indicates resonance of order m and there is an accurate independent determination of these 
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TABLE 18 

NORMALIZED COEFFICIENTS FOR THE GEM-6 MODEL (x 10 6 ) 

ZONALS 


INDEX 

VALUE 

INDEX 

VALUE 

INDEX 

VALUE 

INDEX 

VALUE 

INDEX 

VALUE 

N M 


N 

M 


N 

M 


N 

M 


N 

M 


2 0 

—484.1661 

3 

O 

0.9607 

4 

O 

0.6382 

5 

0 

0.0661 

6 

0 

—0.1461 

7 0 

0.0961 

8 

0 

0.0426 

9 

0 

0.0264 

10 

0 

0.0606 

11 

0 

-0.0528 

12 0 

0.0306 

13 

0 

0.0470 

14 

0 

-0.0206 

15 

0 

-0.0045 

16 

0 

-0.0077 

17 0 

22 0 

0.0192 

-0.0138 

18 

0 

0.0091 

19 

0 

0.0044 

20 

0 

0.0143 

21 

0 

-0.0098 


SECTORIALS AND TESSERALS 


INDEX 

VALUE 

INDEX 

VALUE 


INDEX 

VALUE 

N 

M 

C 


N 

M 

Z 


N 

M 


7 

2 

1 

-0.0009 

-0.0032 

2 

2 

2.4251 

-1.3883 

3 

1 

2.0021 

0.2482 

3 

2 

0.9332 

-0.6311 

3 

3 

0.6969 

1.4260 

4 

1 

-0.5403 

-0.4648 

4 

2 

0.3461 

0.6695 

4 

3 

0.9655 

-0.2073 

4 

4 

-0.1636 

0.3051 

5 

1 

-0.0684 

-0.0842 

5 

2 

0.6651 

-0.3112 

5 

3 

-0.4656 

-0.1947 

5 

4 

-0.2485 

0.0360 

5 

5 

0.1845 

-0.7119 

6 

1 

-0.0734 

- 3 Cl 59 

6 

2 

0.0643 

-0.3740 

6 

3 

0.0115 

0.0098 

6 

4 

-0 0867 

-0.4655 

6 

5 

-0.2747 

-0.5464 

6 

6 

0.0173 

-0.2627 

7 

1 

0.2501 

0.1385 

7 

2 

0.3463 

0.0875 

7 

3 

0.1988 

-0.1844 

7 

4 

-0.2807 

-0.1408 

7 

5 

0.0265 

0.0228 

7 

6 

-0.3074 

0.1213 

7 

7 

0.0624 

0.0048 

8 

1 

0.0102 

0.0579 

8 

2 

0.0610 

0.0860 

8 

3 

-0.0378 

-0.0667 

8 

4 

-0.2311 

0.0284 

8 

5 

—0.0570 

0.0622 

8 

6 

-0.0947 

0.2528 

8 

7 

0.0658 

0.0819 

8 

8 

-0.0832 

0.0701 

9 

1 

0.1426 

0.0137 

9 

2 

0.0552 

-0.0216 

9 

3 

-0.1299 

—0.0727 

9 

4 

-0.0125 

-0.0147 

9 

5 

-0.0036 

-0.0686 

9 

6 

0.0163 

0.1267 

9 

7 

-0.0566 

-0.0037 

9 

8 

0.2519 

-0.0101 

9 

9 

-0.0275 

0.0873 

10 

1 

0.0927 

-0.1343 

10 

2 

-0.0419 

-0.0703 

10 

3 

-0.0383 

-0.0998 

10 

4 

-0.0699 

-0.1247 

10 

5 

-0.0525 

-0.0377 

10 

6 

-0.0550 

-0.1342 

10 

7 

0.0174 

-0.0237 

10 

8 

0.0473 

-0.1213 

10 

9 

0.0960 

-0.0749 

10 

10 

0.1365 

-0.0391 

11 

1 

-0.0087 

0.0431 

11 

2 

-0.0123 

-0.1176 

11 

3 

-0.0334 

-0.0841 

11 

4 

-0.0263 

-0.1033 

11 

5 

0.0842 

0.0406 

11 

6 

-0.0421 

-0.0353 

11 

7 

0.0041 

-0.1123 

11 

8 

—0.0204 

0.0714 

11 

9 

-0.0349 

0.0393 

11 

10 

-0.0464 

-00333 

11 

11 

0.0696 

-0.0356 

12 

1 

-0.0717 

—0.0477 

12 

2 

-0.0530 

0.0621 

12 

3 

0.0G94 

0.0371 

'.2 

4 

-0.0488 

-0.0158 

12 

5 

0.0624 

0.0230 

12 

6 

0.0631 

0.0280 

12 

7 

-0.0261 

0.0127 

12 

8 

-0.0187 

-0.0031 

12 

9 

-0.0002 

0.0251 

12 

10 

0.0231 

-0.0012 

12 

11 

0.0066 

0.0359 

12 

12 

-0.0123 

-0.0103 

13 

1 

-0.0157 

-0.0216 

13 

2 

-0.0454 

-0.0867 

13 

3 

-0.0460 

0.0454 

12 

4 

0.0298 

-0.0670 

13 

5 

0.0586 

0.0469 

13 

6 

-0.0848 

0.0577 

13 

7 

-0.0414 

-0.0482 

13 

8 

-0.0055 

-0.0347 

13 

9 

0.0271 

0.0588 

13 

10 

-0.0240 

-0.0044 

13 

11 

-0.0576 

-0.0830 

13 

12 

-0.0261 

0.0991 

13 

13 

-0.0543 

0.0722 

14 

1 

-0.0038 

0.0480 

14 

2 

-0.0150 

0.0429 

14 

3 

0.0653 

0.0032 

14 

4 

0.0019 

0.0010 

14 

5 

—0.0144 

-0.0216 

14 

6 

0.0166 

-0.0442 

14 

7 

0.0420 

0.0030 

14 

8 

-0.0007 

-0.0414 

14 

9 

0.0140 

0.0552 

14 

10 

-0.0380 

-0.0797 

14 

11 

0.0614 

-0.0313 

14 

12 

0.0047 

-0.0413 

14 

13 

0.0211 

0.0281 

14 

14 

-0.0448 

-0.0016 

15 

1 

0.0333 

-0.0224 

15 

2 

0.0370 

-0 0641 

15 

3 

-0.0457 

0.0279 

15 

4 

-0.0070 

0.0163 

15 

6 

0.0136 

0.0358 

15 

6 

-0.0130 

-0.1076 

15 

7 

0.0751 

0.0668 

15 

8 

-0.0261 

-0.0242 

15 

9 

0.0116 

0.0385 

15 

10 

0.0351 

-0.0480 

15 

11 

-0.0090 

-0.0106 

15 

12 

-0.0338 

0.0145 

15 

13 

-0.0191 

- 0.0000 

15 

14 

0.0036 

-0.0189 

15 

15 

-0.0444 

0.0356 

16 

1 

0.0321 

-0.0091 

16 

2 

-0.0200 

0.0639 

16 

3 

-0.0083 

-0.0205 

16 

4 

0.0252 

0.0306 

16 

5 

0.0120 

0.0173 

16 

6 

0.0321 

0.0136 

16 

7 

-0.0003 

-0.0286 

16 

8 

-0.0456 

-0.0046 

16 

9 

—0.0652 

-0.0676 

16 

10 

-0.0118 

0.0386 

16 

11 

0.0189 

-0.0078 

16 

12 

0.0197 

-0.0195 

16 

13 

0.0034 

-0.0139 

16 

14 

-0.0172 

-0.0430 

16 

15 

-0.0475 

-0.0378 

16 

16 

-0.0376 

-0.0119 

17 

12 

0.0137 

-0.0012 

17 

13 

0.0145 

0.0204 

17 

14 

-0.0111 

-0.0013 

18 

12 

-0.0636 

-0.0269 

18 

13 

-0.0137 

—0.0580 

18 

14 

-0.0127 

-0.0006 

19 

12 

-0.0309 

-0.0300 

19 

13 

-0.0205 

-0 0291 

19 

14 

-O.OOOR 

n rvn i 

20 

12 

-0.0056 

-0.0154 

20 

13 

0.0114 

-0.0282 

20 

14 

0.0055 

-0.0098 

21 

12 

-0.0324 

-0.0175 

21 

13 

-0.0241 

0.0108 

21 

14 

0.0094 

0.0078 

22 

12 

-0.0435 

-0.0065 

22 

13 

-0.0324 

-0.0151 

22 

14 

-0.0077 

-0 0007 
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resonant coefficients, the GEM-6 mode] can be used after replacement of the model 
coefficients by the independent coefficients. 


3. 3.2. 3 12 to 24 Hour Orbital Periods; Low-Order Resonance 

When orbital periods are in the 12- to 24-hour range, the GEM-6 model can be reduced 
considerably in size to accommodate the orbit. The exact size of the model to be used 
depends on the orbital eccentricity and accuracy requirements. For nearly circular 
synchronous (24 hour) orbits, coefficients through (3,3) are usually acceptable. 

Table 19 gives coefficients which may be used to replace corresponding terms in the GEM-6 
model when low-order resonance effects are important. The satellite orbits given in table 20 
were used in computing the values in table 19. 


3. 3.2.4 Polar Orbit Model 

For polar orbits (inclination near 90 degrees), the GEM -6 model can be used. 

3. 3. 2.5 Small Models 

A requirement for a small model with general applicability can be met by truncating the 
GEM-6 model. 


TABLE 19 

LOW-ORDER RESONANT COEFFICIENTS (x 1 0 6 ) 


n 

m 

C 

nm 

S nm 

2 

2 

2.432 

-1.407 

3 

1 

1.878 

0.247 

3 

3 

0.703 

1.470 

4 

2 

0.335 

0.671 

4 

4 

-0.137 

0.374 


1 References 45 and 37. 
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TABLE 20 

24 HOUR SATELLITE OBSERVATIONS USED IN SOLUTION 
OF RESONANT COEFFICIENTS IN TABLE 19 


Longitude Span, 
AX 

Longitude Span, 

AX 

(Degrees East 

Inclination, i 

(Degrees East 

Inclination, i 

of Greenwich) 

(degrees) 

of Greenwich) 

(degrees) 

302-305 

33.0 

210-211 

1.2 

296-301 

32.8 

288-312 

0.6 

196-243 

32.5 

287-288 

0.3 

72-189 

32 ;3 

313-315 

0.3 

65-68 

31.8 

349-352 

1.0 

65-86 

29-32.0 

346-349 

1.0 

178-180 

0.0 

179-194 

1.0 

174-181 

0.0 

324-332 

0.5 

165-172 

0.5 

321-331 

0.8 

160-161 

1.3 

253-256 

2.6 

146-160 

2-3.0 




3.4 Moon 

3.4.1 Mass and Gravitational Constant 

The mass of the Moon is given accurately only in terms of the mass of the Earth: 

M e 

m " 81.302 

When treating the Moon as a point mass, the recommended value of gravitational constant 
of the Moon should be assigned the value: 

M m =4902.78 km 3 /s 2 

In the more detailed lunar gravity models, however, the value of is assigned in each case 
as part of the solution. 

3.4.2 Triaxial Model 

To define the values of the harmonic coefficients of a triaxial Moon, C 20 and C 22 are taken 
as 


C 2Q =-2.07108 X 10' 4 
C 22 =0.20716 X 10' 4 
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3.4.3 Applications 


The triaxial ellipsoid model is recommended for general application except for low altitude 
lunar satellites. For low altitude orb iters, the considerations of section 2. 3. 3. 2 are pertinent. 


3.5 Asteroids and Natural Satellites 

3.5.1 Asteroids 

Although there is little quantitative information about the physical properties of the 
asteroids, orbits and ephemerides for about 1800 of them are known in some detail. 
Ephemerides through the year 2000 have been computed by Duncombe (ref. 59) for the 
four most widely-observed asteroids, Ceres, Pallas, Juno, and Vesta. Reference 3 lists several 
references that can be consulted for detailed ephemerides of the asteroids. 

The masses of the asteroids can be estimated reliably only by observing the perturbed 
motions of other bodies in their vicinity. The mass of Ceres, the largest asteroid, has 
recently been determined from the motion of Pallas to be 6.7 ± 0.4 X 10' 10 of the solar 
mass by Schubart (ref. 60). The mass of Vesta was found by Hertz (ref. 61) to be 1.20 ± .08 
X 10' 1 0 of the solar mass from the motion of the asteroid Arete. 

Diameters of the four most widely observed asteroids have recently been published by 
Dollfus (ref. 57) and Gehrels (ref. 62). These values are derived from astronomical 
observations made by Barnard in 1902. They are presented in table 21. 


TABLE 21 

DIAMETERS OF PROMINENT ASTEROIDS 


Asteroid 

Diameter (km) 

Ceres 

769 ± 20 

Pallas 

490 ± 25 

Juno 

196+ 25 

Vesta 

400 ± 35 


3.5.2 Natural Satellites 

The masses and diameters of the natural satellites are presented in table 22. They are based 
on the estimates of Koslavskaya (ref. 63), Dollfus (ref. 57), De Sitter (ref. 64), and Jeffreys 
(refs. 65 and 66). The masses are presented in units of the mass of the planet about which 
each satellite revolves to avoid introducing the uncertainties of the planetary masses. In 
some cases, a mean value is given with a standard error. In other cases, a range of values is 
given that should be interpreted as the interval over which any value is equally probable. 
When no range is given, the mass or diameter value is so uncertain that it was not possible to 
estimate a range of uniform probability. 
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TABLE 22 

MASSES AND DIAMETERS OF THE NATURAL SATELLITES 




Planet 

Satellite 

Mass 

(Fraction of Planet) 

Diameter (km) 

Earth 

Moon 

0.0123 

3476 ± 2 

Mars 

Phobos (1) 

2.7 x 10 8 

22.46 ± 2 


Deimos (II) 

4.9 x 10 9 

13.4 ± 2.7 

Jupiter 

lo (1) 

(4.696 ±. 06) x 10 5 

3500 ± 150 


Europa (II) 

(2.565 ± .06) x 10 5 

31 10 ± 150 


Gamymede (III) 

(7.845 ± .08) x 10 5 

5550 ± 130 


Callisto (IV) 

(5.603 ±.1 7) x 10 5 

500 ± 150 


V 

2.2 x 10 9 

88 to 160 


VI 

8.5 x 10 10 

64 to 1 84 


VII 

3.5 x 10 11 

24 to 64 


VIII 

7 x 10 13 

6.4 to 18 


IX 

1.5 x 10 12 

7 to 20 


X 

1.0 x 10 12 

7 to 20 


XI 

2.0 x 10 12 

24 


XII 

7.0 x 10' 1 3 

6 to 18 

Saturn 

Mimas (1) 

(6.64 ±0.10) x 10 8 

450 ± 665 


Enceladus (II) 

(1.4 ± 0.47) x 10 7 

550 ± 300 


Tethis (III) 

(1.118 ±0.015) x 10 6 

1 200 ± 200 


Dione (IV) 

(1.90 ±0.07) x 10 6 

820 ± 400 


Rhea (V) 

(3.8 ±3.8) x 10 6 

1300 ±300 


Titan (VI) 

(2.425 ± 0.020) x 10 4 

4850 ± 300 


Hyperion (VII) 

8 x 10 8 

980 to 500 


lapetus (VIII) 

(3.20 ±0.74) x 10' 6 

1 150 ± 100 


Phoebe (IX) 

5 x 10 8 

190 to 540 

Uranus 

Ariel (1) 

28 x 10 6 

760 to 2170 


Umbriel (II) 

8 x 10r® 

500 to 1410 


Titania (III) 

49 x 10 6 

910 to 2600 


Oberon (IV) 

38 x 10 6 

830 to 2380 


Miranda (V) 

1.5 x 10-® 

280 to 820 

Neptune 

Triton (1) 

(1.34 ±0.23) x 10 3 

3770 ± 1500 


Nereid (II) 

3.33 x 10 7 

280 to 800 
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APPENDIX A 

DEFINITION OF SYMBOLS* 


a 

K 

3 e 

A 

b 

B 

C 


e 

f 

F 

g 

G 

i 

I 

K 

K 

K 

m 

n 

n 

P 


Semimajor axis 

Mean equatorial radius of a massive body 
Equatorial radius of a massive body 
Smallest of a body’s principal moments of inertia 
Polar radius of an ellipsoid of revolution (2.2) 

Intermediate of a body’s principal moments of inertia 
Largest of a body’s principal moments of inertia 
Dimensionless spherical harmonic coefficient 
Eccentricity 

Flattening of an ellipsoid of revolution 
Gravitational force 
Gravitational acceleration 

Gravitational acceleration at the equator of a massive body 
Universal gravitational constant (2. 1 ) 

Orbit inclination relative to equator of attracting body 
Moment of inertia relative to any specified axis 
Zonal harmonic of degree n (2.3) 

Factor in computation of K nm (=1 for m=0, =2 for m=£0) (appendix F) 
Normalization factor (appendix F) 

Index connoting the order of a harmonic 
Index connoting the degree of a harmonic 
A unit vector normal to some surface (2.5) 

Some arbitrary point in space 

Legendre polynomial of degree n and order m 

Radial distance from center of mass in spherical coordinates 

Equatorial radius of an equipotential surface near an ellipsoidal body (2.2) 

Dimensionless spherical harmonic coefficient 


♦Numbers in parens give section where symbol used. 
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t Time 

U Gravitational potential due to zonal terms 

V(r, 8, X) Spherical harmonic potential function (Appendix D) 

Y n (0, X) A surface harmonic of degree n (Appendix D) 


Subscripts 

e 

e 

m 

m 

n 

o 

P 

s 


Earth 

Equator 

Moon 

Order (when preceded by n) 

Degree 

Measured at some epoch or reference 

Planet 

Sun 


Greek Letters 

a n Coefficients of degree n used in describing a spheroidal surface (appendix C) 

8 Colatitude (=7t/2 - 0) 

X Longitude 

p Gravitational constant (= GM) 

7r 3.14159265 


p(r) Radius dependent part of the spherical harmonic potential (appendix D) 

0 Latitude 

^nm 5 ^nm Components of the spherical harmonic potential (Appendix D) 
co Rotation rate of a massive body 
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APPENDIX B 

ASTRODYNAMIC CONSTANTS AND UNITS 


The basis of the definitions of the standard astrodynamic constants is Kepler’s laws of 
planetary motion (ref. 67). Analytically stated, Kepler’s Third Law (refs. 67 and 68) says 
that for any planet, the square of its period P is proportional to the cube of its mean 
distance from the Sun, or symbolically, 


where a 



is the semimajor axis of the planet’s orbit about the Sun 


(B-l) 


M s is the Sun’s mass 

M p is the planet’s mass, including planet, atmosphere, and satellites 
k =y/G is the Gaussian gravitational constant 


Historically, equation (B-l) is the basis for the definitions of presently adopted 
astrodynamic units and constants; it contains units of time, mass, and length, and a 
“universal constant” k whose value can be theoretically derived from perfect determinations 
of the observable quantities. Therefore, a convenient set of units was adopted: 


• The unit of mass was taken to be the Sun’s mass, i.e., M s was set equal to unity 

• The unit of time was taken to be a mean solar day, i.e., P was set equal to 
365.2563835* 


• The unit of length was taken to be the mean distance between the Earth and Sun and 
was called the astronomical unit (AU) 

With these units, a value of k could be inferred from the ratio M e /M s where M e is the mass 
of the Earth-Moon system. Reference 67 notes that at the time k was computed, M e /M was 
thought to be (354, 710)' 1 , and k was computed to be 0.01720209895. If the units of 
measurement were to be held fixed, then any refinement in M e /M s would necessitate 
changing k as well. For practical astronomy, however, a consistent value of k is virtually a 
necessity so that k is fixed at the foregoing value by convention. As a result of this 
convention, the unit of length AU equals the mean radius of the Earth’s orbit only if the 
value of M g /M s that was used in deriving k were correct. On the basis of current values of 
M e /M s , the mean radius of the Earth’s orbit a e is 1 .00000003 AU (ref. 3). 

*P actually also includes the time correction 1.1 x 10'^ T where T is measured in centuries since 1 900. 
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APPENDIX C 

GRAVITATIONAL FIELD OF A SPHEROIDAL BODY 
C.1 Mathematical Description 

The derivations, formulas, and notation used here are essentially those of Kaula (ref. 26). 
The expressions for the potential and the radius of an equipotential surface are given in 
terms of the Legendre polynomials P 20 and P 40 (sec. C.2) by 



3 . 2 . 1 

= I sin *’7 



35 . 4 15 . 2 3 

— sin 0 - — sin 0 + — 


where 0 is the planetocentric latitude. The potential at a point located on the body’s surface 
at latitude 0 and distance r from the planet’s center is 


U = 


j£ 

r 


X - J 2 



2 



+ 


12 2 2 
— to r cos 


0 


(C-l) 


where a e is the equatorial radius of the planet 

to is its rotation rate (for an orbiting spacecraft, the rotational 
term is not present in the potential) 

The shape of an equipotential surface for the spheroidal body is symmetric about the axis of 
rotation and is defined by 


r = r 

o 


(1 


+ a 2 P 20 +a 4 


40 


) 


where 


r is the radial distance from the body’s center to a point 
on the equipotential surface with latitude 0 

r Q is the equatorial radius of that surface 
a 2 and a 4 are coefficients to be determined 
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By rewriting r in terms of <j>, and using the binomial expansion for r n , one can write U in 
terms of powers of sin 2 0. The coefficients a 2 and a 4 may be expressed in terms of f and oo, 
whence after combining coefficients of powers of sin 2 0 in the expression for U, the three 
parameters ju, J 2 , and J 4 may be obtained to order f 2 : 


u — a g 
r e e 


1 - f + 


(1 H 

\2 " 4 


m 


T 2 f 1,3 2. 

J 2=3 £ f-J* -J m l 1 -*”” 1 


J 4 = '35 f<7E ‘ 5m> 


where to first order in f, according to Heiskanen and Moritz (p. 74, ref. 69), 


2 

to a 

e 

m 

g e 


and where g e is the gravitational acceleration at the equator. 


C.2. Legendre Polynomials 


Legendre polynomials arise as particular solutions to Legendre’s equation: 

f 2 1 


sin 0 g" (0) + cos 0g' (0) + 


n(n+l) sin 0 - 


sin 0j 


g(0) = 0 


where g(0) is some function of the independent variable 6, n and m are dimensionless 
constants, and ' and " denote first and second derivatives with respect to 8. Making the 
substitutions, t = cos 6 and G(t) = g(0), a solution for G(t) is found to be 
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p 

nm 


(t) 


2 n nl 


(i - t 2 ) m/2 i 


n+m 


dt 


n+m 


(t 2 - l) n , (m/0) 


P 

n 


(t) = P 


no 


(t) 


„n , 1A n 
2 n. dt 


(m = 0) 


where t is cos 6 and 6 is the colatitude. 

The functions P n (t) are polynomials in t. These polynomials may also be obtained from the 
recursion formula 


P„( t )=^P n 2 (t) 

n n n-z 


+ tp (t) 

n n-1 


P o (t) = 1 


by making use of the trigonometric identities: 


2 

cos 


0=jcos 20 +y 


3 

cos 


e 


1 3 

= — cos 36 + — cos 
4 4 


6 , 


etc. 


The Legendre polynomials with m = 0 are summarized in table C-l. 
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TABLE C-1. 

LEGENDRE POLYNOMIALS OF ZERO ORDER 


n P (cos 0) 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


1 

cos 0 

\ (3 cos 2 0-1) 

\ (5 cos 3 0-3 cos 0) 

(35 cos 4 0 - 30 cos 2 6 + 3) /8 
(63 cos 3 0-70 cos 3 0 + 15 cos 0 ) /8 
(231 cos 6 0 - 315 cos 4 0 + 105 cos 2 0-5) /16 
(429 cos 7 0 - 693 cos 5 0 + 315 cos 3 0 - 35 cos 0) /16 
(6435 cos 3 0 - 12012 cos 6 0 + 6930 cos 4 0 - 1260 cos 2 0 + 35)/l28 
(12155 cos 3 0 - 25740 cos 7 0 + 18018 cos 5 0 - 4620 cos 3 0 + 315 cos 0)/l28 
(46189 cos 10 0 - 109395 cos 8 0 + 90090 cos 6 0 - 30030 cos 4 0 + 3465 cos 2 0 - 63)/256 


Associated Legendre functions can be derived from Legendre polynomials by means of the 
equation 


P 

nm 


(t) = (l - t 2 ) m/2 


d m P (t) 
n x ' 


dt 


m 


Associated Legendre functions through n = 10 are given in Table C-2. 
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TABLE C-2. 

ASSOCIATED LEGENDRE FUNCTIONS (1 of 5) 


n 

m 

P (cos 0) 
nm 

i 

i 

sin 0 

2 

i 

3 sin 0 cos 0 

3 

i 

. „ / 15 2 3 \ 

sine( T cos 6-jj 

2 

2 

3 sin 2 0 

3 

2 

2 

15 sin 0 cos 0 

3 

3 

15 sin 2 0 

4 

1 

. /35 3 . 15 

sin 0 cos 0 - — cos 0) 

4 

2 

. 2 „ / 105 2 . 15 \ 

| sin 0 cos 0 - — ) 

4 

3 

3 

105 sin 0 cos 0 

4 

4 

105 sin 4 0 

5 

1 

/ 315 4 . 105 2 15\ 

sin 0 (“g-- cos 6 — ~ cos 0 +- 8/ 

5 

2 

2 / 3 15 3 „ 105 _\ 

sin 0 \— — cos 0 — cos 0 j 

5 

3 

. 3 ^ / 945 2 n 105 \ 

sin 0^— cos e-—J 

5 

4 

4 

945 sin 0 cos 0 

5 

5 

945 sin 5 0 
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TABLE C-2. 

ASSOCIATED LEGENDRE FUNCTIONS (2 of 5) 






TABLE C-2. 

ASSOCIATED LEGENDRE FUNCTIONS (3 of 5) 












TABLE C-2. 

ASSOCIATED LEGENDRE FUNCTIONS (4 of 5) 
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TABLE C-2. 

ASSOCIATED LEGENDRE FUNCTIONS (5 of 5) 


n 

m 

P (cos 0) 
nm 

10 

5 

. 5 „ / 43648605 5 „ 11486475 3 „ 2027025 A \ 

sin 0 1 g cos 0- - cos 0+ cos 0j 

10 

6 

.6^/218243025 4 A 34459425 2 n 2027025\ 

sin 01 g cos 0 - cos 0 + g 1 

10 

7 

. 7 / 218243025 3 n 34459425 .\ 

sin 0 ^ cos 0 cos 0 J 

10 

8 

i 

.8^/654729075 2 A 34459425\ 

sin 0^ 2 cos 0- 2 ) 

10 

9 

654729075 sin 9 0 cos 0 

10 

10 

654729075 sin 10 0 
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APPENDIX D 

SPHERICAL HARMONICS AND ORTHOGONALITY 


The general solution for the gravitational potential exterior to a massive body satisfies 
Laplace’s equation (ref. 67). 


The expression for V in terms of spherical harmonics may be derived from the expression 
for Laplace’s equation in spherical coordinates: 


2 2 Tr 9 / 2 SV 

r V V = — I r — 


Sr \ ar / sin 0 a0 


^-(cos 0 


a v 
a0 


3 2 V 


2 2 
cos 0 a X 


= 0 


(D-l) 


where r = radial distance from center of mass 
0 = colatitude (7t/2 - latitude) 

X = longitude (eastward) 

To solve Laplace’s equation, the variables r, 0, X are separated by the form: 


V(r, 0, X) =p(r) Y(0, X) 


(D-2) 


whence, applying the boundary condition *2^ (V) = 0, the solutions 


p (r) = A r 
n n 


“(n+1) 


and 


(D-3) 


Y (0 > X)-r C' cos mX + s' sin mXJ P (cos 0) 
m=0 \ / 

are found where A , C' , and S' are coefficients whose values are to be determined; 
P nm (cos 6 ) are the Legendre functions (appendix C); and the functions Y n (0, X) are 
known as surface harmonics. The general solution for V is given as a summation of the 
above solutions for P and Y : 

n n 

co A n 

V(r, 0, X) = V — V (C* cos mX + S' sin mX)P (cos 0) (D-4) 

n+1 ' \ nm nm / nm 

n=0 r m=0 
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where V(r, 6 , X) is the gravitational potential in spherical coordinates (sec. 2.3.2 and 
appendix F). Note that in this formulation the product A 0 C Q0 is equivalent to ju as defined 
in section 2. 1 . 

The individual components of the surface spherical harmonics are represented by 


$ C = C cos mX P 
nm nm nm 


(cos 0) 


and 


g 

$ = S sin P 

nm nm nm 


(cos 0) 


An important property of the functions <f> 9 m is that they are orthogonal, that is, 

the integrated product of two different functions, <h° m and over the surface of the 

sphere, a, is zero. 


// a Km < 9 ’ X > < e ' X > dtr = 0 

f/X* 9 ’ X) < r <e ’ x,dii=0 


if s/n or r / m 


(D-5) 


and the integral of the product of two identical functions (s = n and r = m) is 


//.[■ 

//.[• 


(9, X) 
no 


da 


■ //] 


^no < 9 > X > 


* < 0 ’ X) 
nm 


da= //[ 


* m (6, X) 
nm 


2 da = -fr 

n + 1 


2 2tt (n + m)l 

a 2n + 1 (n - m)l 


(D-6) 


Orthogonality makes the spherical harmonics the natural means of representing a function 
over a spherical surface (analogous to the use of Fourier series for functions in a rectilinear 
space). The orthogonality property means that the effect of each term of the harmonic 
series is unique. Thus terms or groups can be studied as independent subsets of the complete 
series. Further, the spherical harmonic expansion model is not limited to the Earth; it has 
been applied to the Moon and some of the planets. 

In the case of the Earth, the spherical harmonic model is discussed in section 3.3. 
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APPENDIX E 

SPHERICAL HARMONIC REPRESENTATION OF 
THE GRAVITATIONAL POTENTIAL 


E.l. Zonal Harmonics 

The expression below: 

(C nm cosmX + S nm sinmX ) P nm < COS 

is a surface spherical harmonic of degree n and order m. If m = 0, the harmonic is a constant 
multiple of the Legendre function P no (cos 0). The function, P no (cos 0), has n distinct 
zeros between 0 =0 and 0 = tt (between -ir/2 and it/ 2 latitude) arranged symmetrically about 
0 = 2 * (% E-l). 



Figure E-1a.- P no (cos 0) for n Even 
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In a coordinate system having its origin at the center of a sphere, the function P n0 (cos 0 ) 
vanishes on n circles of constant latitude (fig. E-2). Similarly, the locus of points on the 
sphere at which the function has a constant value consists of a number of parallel circles. 
Because of this division of the sphere into latitude zones in which P no (cos 8 ) alternates in 
sign, the functions, P no (cos 0), are called zonal harmonics. 

Many authors have analyzed the zonal harmonics of the Earth’s gravitational potential, e.g., 
Kozai (ref. 70) and King-Hele (ref. 71) have represented the potential at a point, (r, 0), by 



-E J 

n=0 


/ a \ n 

' ~ 

n\ r / 


(cos 0) 


(E-l) 


In what follows, the standard convention of omitting the o subscript from P no will be 
adopted, so that, in general J n = -C . 

In pre-1958 literature (refs. 24 and 25), only even order zonal harmonics were believed to 
be significant in the shape and potential of the Earth. This assumption would be true if the 
Earth were an equilibrium figure of rotation. Even though it was admitted that this 
assummption was not perfect, the Earth was believed to be sufficiently near equilibrium to 
make all odd J n negligible. Henriksen (ref. 72) and Cook (ref. 36) were among the early 
researches in the interpretation of satellite results to accept the possibility that the 
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Figure E-2. — Alternate Positive and Negative Regions of P 40 (cos d) = 0 
(Shaded Areas are Positive) 


flattening of the Earth is not necessarily that of a body in perfect pressure equilibrium; this 
led to the adjustment of the value of the second harmonic. O’Keefe, Eckels, and Squires 
(ref. 73) were the first to report the existence of odd zonal harmonics. They attributed the 
80-day periodic variation of the eccentricity in the orbit of Vanguard 1 to the presence of 
the third zonal harmonic in the Earth’s field. The announcement of this result gave rise to 
the term “pear-shaped” in describing the Earth’s shape. 

Zonal harmonics of even degree give rise to secular perturbations of the orbital elements 12 
and co, and both even and odd zonals give rise to long period perturbations of e, i, 12 and to. 
Therefore, their influence can be detected in changes of orbital parameters that are 
integrated over many revolutions of satellites. The discussion in section 2.2 of the 
relationship between the zonals J 2 and J 4 and the secular changes in 12 and to may be 
generalized to include numbered zonals of higher degree (ref. 27). 

E.2 Tesseral Harmonics 

If 0 < m < n, then the associated Legendre functions, P nm (cos 6), change their sign n-m 
times in the colatitude interval 0 < 6 < 7r and the surface harmonic representation takes the 
form 

fc cos mX + S sin m\) P (cos 0) 

\ nm nm / nm 

The functions cos mX and sin mX have 2m zeros in the longitude interval 0 < X < 27r. The 
geometrical representation of such a harmonic is shown in figure E-3 where the sphere is 
divided into compartments (tesserae) which are alternately positive and negative; the 
harmonics are called tesseral harmonics. 
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Figure E-3. - Alternate Positive and Negative Regions of 
P 10 5 (cos 0) {C 10 5 cos 5\ + S 10 5 sin 5X) 
(Shaded Areas are Positive) 


The perturbations caused by tesseral harmonics usually have short periods (sub-multiples of 
the planet’s period of rotation). 

E.3 Sectorial Harmonics 

If m = n, then the surface harmonic takes the form 


(c cos n\ + S sin n\ ) sin 11 0 
\ nn nn ' 

which is represented on a sphere in figure E-4. In this form the harmonic oscillates in sign 
within longitude bands separated by n great circles passing through 8 = 0 and 8 = n, i.e., the 
poles. Because the sphere is thus divided into 2n sectors, the designation, sectorial 
harmonics, is used for this type of harmonics. Except for resonance effects, perturbations 
on satellite motion that result from the sectorial harmonics have short periods (submultiples 
of a day). The determination of the sectorial harmonics is similar to that of tesseral 
harmonics. 
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Figure E-4. — Alternate Positive and Negative Sectors for 
(C 55 cos 5X + S 55 sin 5A) sin 50 
(Shaded Areas are Positive) 


E.4 Resonant Harmonics 

The orbits most nearly commensurate with the Earth’s rotation have 12- and 24-hour 
periods. The constraints placed on the low-order harmonics of the geopotential by analysis 
of these orbits were presented by Wagner in reference 74. The high altitude resonant data 
provide verification of the low-order gravity terms which reflect average dynamic effects 
over large areas. In the same way, surface gravity data over limited regions give strength to 
the terms of high order and degree (ref. 26). 

The actual number of orbits significantly perturbed by resonance phenomena is much 
greater than formerly supposed; this is an important consideration for mission planners 
concerned with tracking of Earth satellites. The altitudes and periods of Earth resonant 
orbits are listed in table E-l. Likewise, resonance may be important for planning orbits of 
other planets. 

The discovery of the resonance of satellite orbits with the gravitational potentials has 
provided a means of obtaining values for harmonic coefficients whose contributions, 
otherwise, might have been too small to detect. The main investigators in the field have been 
Wagner, working mainly with low-order resonant harmonics; Gaposchkin (ref. 43), working 
with 9th, 12th, 13th, and 14th order; Anderle (ref. 40), with 13th order; Yionoulis (ref. 
75), 13th order; Douglas and Marsh (ref. 76), with 13th order; and King-Hele, et al. (ref. 77) 
with 1 5th order resonant coefficients. 


E.5 The Gravitational Potential 
Using Spherical Harmonics 


The general expression for the gravitational potential at a point (r, 6 , X), external to a planet, 
including zonal, tesseral, and sectorial harmonics is 



cos m\ + S 

nm 


sin mX 


P (cos 0) 
nm 


(E-2) 


where a e is the equatorial radius of the planet. 

Equations D-4 and E-2 are different because in equation D-4 the integration constants have 
been written in terms of a e and p (the central force coefficient) and in equation E-2 the 
harmonic coefficients C and S„ m are dimensionless. In equation E-2, terms of degree 1 
are omitted because they represent terms that result in an offset of the center of mass of the 
body relative to the center of the coordinate system to which the expression for U is 
referenced. In practice, such offsets are usually negligible or zero. In addition, the terms C 2 j 
and S-, j are set to zero because of the coincidence of the Earth’s rotation axis and its axis of 
maximum moment of inertia. 


TABLE E-1 

ALTITUDES AND PERIODS OF RESONANT SATELLITE ORBITS 


Resonant 

Order 

m 

Satellite Altitude 

Period (min) 

Naut. Miles 

km 

1 

1436.0 

19,320 

35,781 

2 

718.0 

10,895 

20,178 

5 

287.2 

3,340 

6,186 

8 

179.5 

2,248 

4,163 

9 

159.56 

1,816 

3,363 

10 

143.6 

1,461 

2,706 

11 

130.05 

1,146 

2,122 

12 

119.67 

899 

1,665 

13 

110.5 

671 

1,243 

14 

102.57 

475 

880 


56 



APPENDIX F 

NORMALIZATION AND CONVERSION 


F.1 Normalization 

The conventional harmonics C nm and S nm are sometimes replaced by normalized 
harmonics, C nm and S nm . Normalization has the effect of giving the coefficients C nm and 
S a clear physical interpretation. For any two harmonic perturbations that are observed 
to have the same magnitude, the following expression has the same value 




c 2 


nm 


+ S 


2 


nm 


This relationship is seen more clearly in the formulation of the potential in terms of the 
coefficients J™ and X™ where the magnitude of the perturbation is completely contained in 
J™ and the phase angle is X™. Thus, 


J m = / c 2 +s 2 ~ 
n a/ nm nm 


and perturbations of equal magnitude would be represented by equal values of J™ . 

The relation between normalized and unnormalized coefficients in the C nm and S nm 
formulation is given in reference 27 by 


C 

nm 


C 

nm 

K 

nm 


S 

nm 


S 


nm 

K 

nm 
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where 


K 

nm 


K(2n + 1) 


(n - m)l 
(n + m)'. 


1/2 


and where 

K = 1 when m = 0 
K = 2 when m / 0 


The selection of this particular expression for K nm may be seen to arise from the results of 
the integrals in equation D-6. A table coefficients for normalization of harmonic coefficients 
is given in table F-l. 


F.2 Conversion 


Committee 7 of the International Astronomical Union (ref. 78) recommended the following 
form and notation for the general expression for the Earth’s gravitational potential: 




(sin 0) 


C cos 
nm 


mX 


+ S sin mX I 
nm 


(F-l) 


where 


M = GM = product of universal gravitational constant times the mass of 
attracting body 

a e = mean equatorial radius of Earth 
r = distance from geocenter to point of observation „ 

0 = geocentric latitude 
X = east longitude 

n,m = indices indicating degree and order 
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labie F-1 

NORMALIZATIONS FACTORS (1 OF 2) 


n m Kfim 


0 

0 

1.000000 

X 

10 

0 

1 

0 

1. 732049 

X 

10 

0 

2 

0 

2.236065 

X 

le 

0 

3 

0 

2. 645748 

X 

le 

0 

4 

0 

2-999997 

X 

IP 

0 

5 

0 

3. 31 6620 

X 

10 

6 

6 

0 

3. 605545 

X 

10 

0 

7 

0 

3. 872976 

X 

10 

0 

e 

0 

4. 123096 

X 

10 

0 

9 

0 

4. 358893 

X 

10 

e 

10 

0 

4. 582563 

X 

10 

0 

11 

0 

4.795625 

X 

10 

0 

12 

0 

4.999991 

X 

10 

0 

13 

0 

5. 196136 

X 

10 

0 

14 

0 

5. 36 5150 

X 

10 

0 

15 

0 

5. 567744 

X 

10 

0 

16 

0 

5. 744547 

X 

10 

0 

17 

0 

5.916055 

X 

10 

6 

18 

0 

6.062737 

X 

10 

0 

19 

0 

6. 244961 

X 

10 

0 

20 

0 

6. 403098 

X 

10 

0 

21 

0 

6. 557402 

X 

10 

0 

22 

0 

6.708173 

X 

10 

0 

23 

0 

6. 855634 

X 

IP 

P 

24 

0 

6. 999969 

X 

16 

0 

25 

0 

7. 141410 

X 

10 

0 

26 

0 

7. 262082 

X 

10 

e 

27 

0 

7. 416167 

X 

10 

0 

28 

0 

7. 549786 

X 

10 

0 

29 

0 

7. 681094 

X 

10 

0 

30 

0 

7.810185 

X 

10 

0 

1 

1 

1. 732049 

X 

10 

e 

2 

1 

1.290996 

X 

10 

0 

3 

1 

1.080125 

X 

10 

0 

4 

1 

0.948683 

X 

10 

0 

5 


0.856350 

X 

10 

0 

6 

1 

0.786795 

X 

10 

0 

7 

1 

0. 731923 

X 

10 

0 

8 

1 

0. 667 186 

X 

le 

0 

9 

1 

0. 649788 

X 

10 

0 

10 

1 

0. 6179P9 

X 

10 

0 

11 

1 

0. 59P321 

X 

1C 

0 

12 

1 

0. 566140 

X 

10 

0 

13 

1 

0- 544706 

X 

10 

0 
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C , S = dimensionless numerical coefficients 

n m n m 

PJJ (sin0) = associated Legendre function 

Variations in form normally occur in the dimensionless numerical coefficients C nm and 
S ; occasionally they are found in the form of V. For example, P™ (cos 9 ) is often used in 
place of P™ (sin 0) where 6 represents colatitude. It may also be noted that the notations 
P m , P , and P are all equivalent. 

The following list gives the potential forms used by various investigators and the required 
scaling and variation for reducing them to the international form (foregoing equation F-l). 


1. Form Used by Moritz, Cambridge Research Laboratory, and Others 




n 


J P 
n n 


(sm 0) 


n 


* E, < J , 

m=l 


cos mX 


+ K sin mX) P m (sin 0) 
nm n 


(F-2) i 


The important difference between this equation and equation F-l is the sign 
ahead of the summation: 


-'I 

- j 

S 

= - K 

n, m 

n, m 

n, m ] 

C 

- J 

s 

= 0 

n,o 

n 

n,o 


n, m 


2. Form Used by Jeffreys, O’Keefe, and Others 



m 


(A cos mX + B sin mX) P 
nm nm n 


(sin 0) 
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n, m 


a a 

n, m \ e i 


3. Form Used by Kozai and Newton and Others 


t i+ z £ (-7) c p > in « 

L n=l m=0 


cos m f X - X 


The conversion equations are 


„ mi . m 

C = J cos m\ 
n, m n n 


_ T m . . m 

S = J sin mX 
n, m n n 


4. Form Used by Mueller and Others 


00 n / a \n+l 

U=T"' V — — ) fa cosmX+b sinmXjP (sin 0) 

L-J n \ r / V nm nm ' n 

M A — A ' 


n=0 m=0 


The conversion is 


C =1— a 
n, m \ (i / n, m 


S =1— b 
n,m \ n I n, m 
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5. Special Form for Zonais Only 

This form is used only when zonal terms are considered adequate to describe the 
potential. The form recommended for this case by committee 7 of the International 
Astronomical Union (Hagihara, ref. 78) is 


U = 




P (sin 0) 


This is the same as the equation F-l without the tesseral and sectorial terms. 


C = - J and 

n,o n 

o 

II 

o 

a 

CO 

Early Form Used by JPL and Others 


„ /a \2 „/ a \ 3 

U'fI 1 — J (-f) P 2 < Sln ®> - s(-T-) p 3 

+ w d {t) V siII « 


The conversion relations are 


C 2,0 = ~3 J 

s„ = 0 
2, 0 

°3,0=-T H 

S = 0 

3,0 

C 4,0 = f D 

O 

II 

o 

CO 


The coefficients in front of J, H, and D are actually combined with the Legendre 
polynomial, which is not explicitly given in the expression for U . 
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7. Form Used by Sterne, Baker, Herrick, and Others 


U = F p 2 <»*“« -fni-r) 3 p 3 < sta « 


+ h k (t) U (sin 0) 


The conversion is 


C 2,0 = -I J 


S 2,0 ° 


C 3,0 " 5 H 


S„ = 0 
3,0 


C 4,0 15 K 


S, „ = 0 

4,0 


Here again the coefficients in front of J, H and K are combined with the Legendre 
polynomial. 

8. Form Used by RCA at Air Force Eastern Test Range 


u 1 " p 2 (Sin 0) + ~^4 p 4 (Sin 0) 
3r 5r 


The conversion relations are 


C 2,0 = ‘ a/ l 3a e 


1 , 0 - 4 %) 


S 2,0 =0 


S 4,0=° 
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9. Form Used by JHU Applied Physics Laboratory 



1 + 


El 

o —o 


ixi r 

P (sin0) C cos m\ + S sin mX 

i [_ 


APL’s unnormalized coefficients are identical to those used in equation F-l where the 
degree index n has been replaced by K. The normalization factors used by APL, however, 
differ from those presented in this appendix. 




APPENDIX G 

PROCEDURE FOR ESTIMATING 
PERTURBATION MAGNITUDES AND TRUNCATING 


The selection of the size of an Earth gravity model is influenced by the magnitudes of the 
perturbations caused by each term. For gravity harmonics of degree n, these will vary with 
r -(n + 1 ) w j iere r j s the geocentric radius. The relative magnitudes of the perturbations caused 
by terms of degree n are shown graphically in figure G-l , which is an expansion of a similar 
figure in reference 79. For purposes of comparison, the gravitational effects of the Sun and 
Moon and the effects of atmospheric drag also are shown. In using the figure, it should be 
assumed that all geocentric radius vectors point in the direction of the Sun and that the 
Moon lies between the Earth and the Sun. The curves of drag perturbations assume circular 
orbital velocity at each altitude. 

As long as the perturbations caused by the individual terms of the spherical harmonic series 
remain small, i.e., conditions of perfect or “deep” resonance are avoided, it is possible to 
obtain expressions. for the perturbations of the orbital elements caused by each V nmpq . 
Assuming that co, M and are all independent of time, Kaula (ref. 26) obtained: 


Aa = 2F G (n - 2p + q) D 
nmpq nmp npq nmpq 


Ae = F G 
nmpw nmp npq 


Ai = F G 
nmpq nmp npq 


where p is the inclination function subscript 

(P = 0,1,2 n). 

q is the eccentricity function subscript 
(q = n - 2p). 

Equation G-l can be used to evaluate the along track, cross track, and normal orbit element 
perturbations caused by each term of the harmonic gravity potential. A computer program 
called HAP is available from the Geodynamics Branch, NASA Goddard Space Flight Center. 
On the basis of equation G-l, HAP provides rapid perturbation estimates for any desired 
orbit of low or moderate eccentricity. 


f 1 - e2 " A - e“ (n - 2p q) - (n - 2p) 


D 


nmpq 

ae 


(G-l) 


(n - 2p) cos (i) - m 


D 


nmpq 


. \fl - sin i 
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NASA SPACE VEHICLE 
DESIGN CRITERIA MONOGRAPHS 


ENVIRONMENT 
SP-8005 
SP-8010 
SP-801 r 
SP-8013 

/ 

SP-801 7 
SP-8020 
SP-8021 
SP-8023 
SP-8037 

SP-8038 

SP-8049 

SP-8067 

SP-8069 

SP-8084 

SP-8085 

SP-8091 

SP-8092 

SP-8103 


Solar Electromagnetic Radiation, revised May 1971 

Models of Mars’ Atmosphere (1974) revised December 1974 

Models of Venus Atmosphere (1972), revised September 1972 

Meteoroid Environment Model— 1969 (Near Earth to Lunar Surface), 
March 1969 

Magnetic Fields— Earth and Extraterrestrial, March 1969 
Mars Surface Models (1968), May 1969 

Models of Earth’s Atmosphere (90 to 2500 km), revised March 1973 
Lunar Surface Models, May 1969 

Assessment and Control of Spacecraft Magnetic Fields, September 
1970 

Meteoroid Environment Model— 1970 (Interplanetary and Planetary), 
October 1970 

The Earth’s Ionosphere, March 1 97 1 

Earth Albedo and Emitted Radiation, July 1971 

The Planet Jupiter (1970), December 1971 

Surface Atmospheric Extremes (Launch and Transportation Areas), 
revised June 1974 

The Planet Mercury (1971), March 1 972 
The Planet Saturn (1970), June 1972 

Assessment and Control of Spacecraft Electromagnetic Interference, 
June 1972 

The Planets Uranus, Neptune, and Pluto (1971), November 1972 
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SP-8105 

SP-8111 

SP-8116 

SP-8117 

SP-8118 

STRUCTURES 

SP-9011 

SP-8002 

SP-8003 

SP-8004 

SP-8006 

SP-8007 

SP-8008 

SP-8009 

SP-8012 

SP-8014 

SP-8019 

SP-8022 

SP-8029 

SP-8031 

SP-8032 

SP-8035 

SP-8040 

SP-8042 

SP8043 


Spacecraft Thermal Control, May 1973 

Assessment and Control of Electrostatic Charges, May 1974 
The Earth’s Trapped Radiation Belts, March 1975 
Gravity Fields of the Solar System, April 1975 
Interplanetary Charged Particle Models (1974), March 1975 

Buffeting During Atmospheric Ascent, revised November 1970 

Flight- Loads Measurements During Launch and Exit, revised June 
1972 

Flutter, Buzz, and Divergence, July 1964 
Panel Flutter, revised June 1972 

Local Steady Aerodynamic Loads During Launch and Exit, May 
1965 

Buckling of Thin-Walled Circular Cylinders, revised August 1968 

Prelaunch Ground Wind Loads, November 1965 

Propellant Slosh Loads, August 1968 

Natural Vibration Modal Analysis, September 1968 

Entry Thermal Protection, August 1968 

Buckling of Thin-Walled Truncated Cones, September 1968 

Staging Loads, February 1969 

Aerodynamic and Rocket-Exhaust Heating During Launch and 
Ascent, May 1969 

Slosh Suppression, May 1969 

Buckling of Thin-Walled Doubly Curved Shells, August 1969 
Wind Loads During Ascent, June 1970 
Fracture Control of Metallic Pressure Vessels, May 1970 
Meteoroid Damage Assessment, May 1 970 
Design— Development testing, May 1970 
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SP-8044 

SP-8045 

SP-8046 

SP-8050 

SP-8053 

SP-8054 

SP-8055 

SP-8056 

SP-8057 

SP-8060 

SP-8061 

SP-8062 

SP-8063 

SP-8066 

SP-8068 

SP-8072 

SP-8077 

SP-8079 

SP-8082 

SP-8083 

SP-8095 

SP-8099 


Qualification testing, May 1970 
Acceptance testing, April 1 970 

Landing Impact Attenuation for Non-Surface-Planning Landers, 
April 1970 

Structural Vibration Prediction, June 1970 

Nuclear and Space Radiation Effects on Materials, June 1970 

Space Radiation Protection, June 1970 

Prevention of Coupled Structure-Propulsion Instability (Pogo), 
October 1970 

Flight Separation Mechanisms, October 1 970 

Structural Design Criteria Applicable to a Space Shuttle, revised 
March 1972 

Compartment Venting, November 1 970 

Interaction with Umblicals and Launch Stand, August 1970 

Entry Gasdynamic Heating, January 1971 

Lubrication, Friction, and Wear, June 1971 

Deployable Aerodynamic Deceleration Systems, June 1971 

Buckling Strength of Structural Plates, June 1971 

Acoustic Loads Generated by the Propulsion System, June 1971 

Transportation and Handling Loads, September 1971 

Structural Interaction with Control Systems, November 1971 

Stress-Corrosion Cracking in Metals, August 197 1 

Discontinuity in Metallic Pressure Vessels, November 1971 

Preliminary Criteria for the Fracture Control of Space Shuttle Struc- 
tures, June 1971 

Combining Ascent Loads, May 1972 
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GUIDANCE AND CONTROL 


SP-8015 

SP-8016 

SP-8018 

SP-8024 

SP-8026 

SP— 8027 

SP-8028 

SP-8033 

SP-8034 

SP-8036 

SP-8047 

SP-8058 

SP-8059 

SP-8065 

SP-8070 

SP-8071 

SP-8074 

SP-8078 

SP-8086 

SP-8096 

SP-8098 


Guidance and Navigation for Entry Vehicles, November 1968 

Effects of Structural Flexibility on Spacecraft Control Systems, 
April 1969 

Spacecraft Magnetic Torques, March 1969 

Spacecraft Gravitational Torques, May 1969 

Spacecraft Star Trackers, July 1970 

Spacecraft Radiation Torques, October 1969 

Entry Vehicle Control, November 1969 

Spacecraft Earth Horizon Sensors, December 1969 

Spacecraft Mass Expulsion Torques, December 1969 

Effects of Structural Flexibility on Launch Vehicle Control Systems, 
February 1970 

Spacecraft Sun Sensors, June 1970 

Spacecraft Aerodynamic Torques, January 1971 

Spacecraft Attitude Control During Thrusting Maneuvers, February 
1971 

Tubular Spacecraft Booms (Extendible, Reel Stored), February 1971 

Spacebome Digital Computer Systems, March 1971 

Passive Gravity-Gradient Libration Dampers, February 1971 

Spacecraft Solar Cell Arrays, May 1971 

Spacebome Electronic Imaging Systems, June 1971 

Space Vehicle Displays Design Criteria, March 1972 

Space Vehicle Gyroscope Sensor Applications, October 1972 

Effects of Structural Flexibility on Entry Vehicle Control Systems, 
June 1972 
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SP-8102 Space Vehicle Accelerometer Applications, December 1972 

CHEMICAL PROPULSION 

SP-8025 Solid Rocket Motor Metal Cases, April 1970 

SP-8039 Solid Rocket Motor Performance Analysis and Prediction, May 197 1 

SP-8041 Captive-Fired Testing of Solid Rocket Motors, March 1971 

SP-8048 Liquid Rocket Engine Turbopump Bearings, March 1971 

SP-805 1 Solid Rocket Motor Igniters, March 197 1 

SP-8052 Liquid Rocket Engine Turbopump Inducers, May 197 1 

SP-8064 Solid Propellant Selection and Characterization, June 1 97 1 

SP-8073 Solid Propellant Grain Structural Integrity Analysis, June 1973 

SP-8075 Solid Propellant Processing Factors in Rocket Motor Design, October 197 1 

SP-8076 Solid Propellant Grain Design and Internal Ballistics, March 1972 

SP-8080 Liquid Rocket Pressure Regulators, Relief Valves, Check Valves, Burst 

Disks, and Explosive Valves, March 1973 

SP-808 1 Liquid Propellant Gas Generators, March 1972 

SP-8087 Liquid Rocket Engine Fluid-Cooled Combustion Chambers, April 1972 

SP-8088 Liquid Rocket Metal Tanks and Tank Components, May 1974 

SP-8090 Liquid Rocket Actuators and Operators, May 1973 

SP-8094 Liquid Rocket Valve Components, August 1973 

SP-8097 Liquid Rocket Valve Assemblies, November 1973 

SP-8 100 Liquid Rocket Engine Turbopump Gears, March 1974 

SP-8 101 Liquid Rocket Engine Turbopump Shafts and Couplings, September 1 972 

SP-81 10 Liquid Rocket Engine Turbines, January 1974 


NASA- Langley, 1975 
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